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Abstract 

The equations of motion for electromechanical systems are traced back 
to the fundamental Lagrangian of particles and electromagnetic fields, via 
the Darwin Lagrangian. When dissipative forces can be neglected the 
systems are conservative and one can study them in a Hamiltonian for- 
malism. The central concepts of generalized capacitance and inductance 
coefficients are introduced and explained. The problem of gauge inde- 
pendence of self-inductance is considered. Our main interest is in magne- 
tomechanics, i.e. the study of systems where there is exchange between 
mechanical and magnetic energy. This throws light on the concept of mag- 
netic energy, which according to the literature has confusing and peculiar 
properties. We apply the theory to a few simple examples: the extension 
of a circular current loop, the force between parallel wires, interacting 
circular current loops, and the rail gun. These show that the Hamilto- 
nian, phase space, form of magnetic energy has the usual property that 
an equilibrium configuration corresponds to an energy minimum. 
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1 Introduction 



Electromagnetism is usually taught at the undergraduate level without men- 
tion of Lagrangians, Hamiltonians, or the principle of least action. In modern 
theoretical physics of gauge field theory, however, the concept of an invariant 
Lagrangian density has become the standard starting point. The Lagrangian 
formalism of analytical mechanics was introduced into electromagnetism already 
by Maxwell in his Treatise T] who, using this approach, derives equations for 
electric circuits and for electromechanical systems. Since then its importance 
has kept growing. One can therefore argue that this set of tools should be better 
known and become accessible at an earlier stage in the physics curricula. This 
review attempts to be an aid in such efforts. 

Our starting point is the basic Lagrangian density of classical electrodynam- 
ics as set down early in the last century by Larmor and Schwarzschild. From 
there we proceed to neglect radiation which leads us to the Darwin Lagrangian 
[3. Then the path to classical linear circuit theory is traced. It is pointed out 
that the similarity between the Lagrangian formulation of mechanics and of cir- 
cuit theory has deep physical reasons and is not just a formal similarity. In a 
long Appendix the generalized capacitance coefficients and the coefficients of self 
and mutual induction of circuit theory are derived, investigated and explained. 

Lagrangians of electromechanical systems are also seen to arise from the 
Darwin Lagrangian by introducing suitable constraints, or assumptions, on the 
possible movements of both the charged particles and the neutral matter in 
the system. We concentrate on magnetomechanical problems, i.e. problems 
where there is a magnetic interaction energy involving macroscopic matter. As 
examples of such problems we consider the extension of a circular loop of current, 
the attraction of parallel currents, the interaction between two circular loops of 
current, and the rail gun. Finally we discuss the properties of the concept of 
magnetic energy and clarify some tricky points. 



2 Lagrangian electrodynamics 

In modern physics one has found that the most reliable and fundamental starting 
point in theoretical investigations is the principle of least action. The action 
is a scalar quantity constructed from a Lagrange density which is a function of 
the relevant particle and field variables and their (normally first) derivatives. 
The action for classical electrodynamics is the time integral of the Lagrangian 
C which has three parts, 

/: = An + A + A. (1) 

The first part is the Lagrangian for free non-interacting particles, 

N N 

= ^ £„, a = 51 -^-^2 (2) 
a— 1 a— 1 

In the non-relativistic approximation, which we mostly assume valid, it is simply 
the kinetic energy. The second is the the interaction Lagrangian, 

A = / (-J -A-gA dV. (3) 
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It was published in 1903 by Karl Schwarzschild (1873 - 1916) and describes the 
interaction of the charge and current density of the particles with the electro- 
magnetic potentials. The third and final part is the field Lagrangian, 

C,^^ j{E^~B^)dV, (4) 

originally suggested by Joseph Larmor (1857 - 1942) in 1900. The connection 
with ^ is via the identifications, 

1 dA 

E = -V(t) — , and B = V X A. (5) 

c ot 

Maxwell's homogeneous equations are identities obtained by taking the curl 
of the first of the equations (O, and the divergence of the second. Maxwell's 
remaining, inhomogeneous equations, and the equations of motion for the par- 
ticles under the Lorentz force, are all obtained from the variation of the action, 
S = J £dt, with £ from ([T]). It is this joining of both the equations determining 
the fields from the sources, and the equations of motion of the sources due to the 
fields, into a single formalism, that is the strength and beauty of this approach. 

The variational approach to electromagnetism outlined above can be found 
in many of the more advanced textbooks on electrodynamics [31 SI O [3 [S] . 
More specialized works are Yourgrau and Mandlestam 'W , Doughty [TU], and 
Kosyakov [TT] . 



2.1 The Darwin Lagrangian 

In many types of problems one can neglect the radiation of electromagnetic 
waves from the system under study, since this phenomenon is proportional to 
c~^. In those circumstances the field Lagrangian C[ can be rewritten and one 
finds that £f = —\£i- Inserting this in H]) we get, 

/: = /:m + \£:, (6) 

for the relevant Lagrangian in the non-radiative case. When the motion of a 
charged particle is known one can find the potentials, 0, A, that it produces, the 
so called retarded, or Lienard-Wiechert potentials. Expanding these to order 
(w/c)^, one finds that acceleration vanishes from the Lagrangian (since it only 
contributes a total time derivative to this order). The result is a Lagrangian 
that contains only particle positions and velocities. There are then no indepen- 
dent electromagnetic field degrees-of-freedom. Everything is determined by the 
positions and velocities of the charged particles, and the resulting Lagrangian is 
the Darwin Lagrangian [2], as derived by Charles Galton Darwin (1887 - 1962), 
a grandson of the great naturalist, in 1920. 
The Darwin Lagrangian can be written, 

£D=£n. + l f (-j -A-qA dV, (7) 



2./ Vc 

i.e. Eq. ([5]), where. 



ir> = / 4^, (8) 



r — r' 
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which is exact in the Coulomb gauge, and where, 

^'2cJ \r-r'\ ^ ' 

Here e^v = {r—r')/\r—r'\. This specific form of tlie Darwin vector potential can 
be traced back to a fairly large retardation effect in the Lorenz gauge Coulomb 
potential. Its effect is included in the Darwin approximation which, however, 
uses a Coulomb gauge, V • A = 0. 

A more familiar form of the Darwin Lagrangian, for N point particles, is 
obtained by introducing, 

N N 

g{r) = ^ ea S{r - ra{t)) , and j{r) = ^ eaVa{t) S{r - ra{t)) , (10) 

a— 1 a—l 

in the expressions ([7]) - (O given above. After skipping self interactions one 
obtains, 

1^/1 \ 

>Cd = £m + 2 X! ( ■ ^a{ra) - ea4'a{ra) \ , (11) 

a—l ^ ^ 

where. 



N 



6(#a) I "I 

and, 



(12) 



N 



^^(r) = ^ ebbb + (^b ■ er,r)e,.,r] ^^^^ 



2c\r ~ rfcl 



Here and t^a are particle position and velocity vectors respectively, ma and Ca 
their rest masses and charges respectively, while e^^r = [i" — '"b)/!^ ^ ''"fcl- There 
are no independent field degrees-of-freedom and hence no gauge invariance in the 
Darwin formalism, which entails action-at-a-distance. Retardation is included 
to order [v/cf' , a fact which is often missed in the literature. 

The Darwin approach to electromagnetism is only briefly mentioned in some 
advanced textbooks [H HI [S] . A book by Podolsky and Kunz ^12^ is a bit more 
thorough. Several good fundamental and pedagogical studies can, however, be 
found in the literature [HIIllIIlllinillTlIIHlIinilin!- 

2.2 The kinetic energy of currents 

In the free particle Lagrangian Cm of Eq. the approximation, 

Cma~-mac^ + '^vl+'^vl (14) 

is usually done, because of the validity of the Darwin approach to order (u/c)^. 
Here we will be concerned with systems in which there are macroscopic charge 
and current densities confined to electrically conducting matter. In 1936 Darwin 
[21j found that the magnetic energy contribution to the inertia of the conduction 
electrons is roughly 10* greater that the contribution from their rest mass. 
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This means that, so caUed, inductive inertia dominates. For the dynamics of 
macroscopic currents and charges in fixed conductors (electric circuit theory) 
one can consequently also neglect the free particle Lagrangian /!ni- In plasma 
physics the neglect of particle inertia is called the force free approximation [22] . 
Skipping Crn, 

£lc - ^ y (^^3 -A-Qcl^ dV, (15) 

is all that then remains of ([7|). This Lagrangian, together with A and (j) given 
by ([8]) and ^ respectively, describes electromagnetic systems with inductive 
and capacitive phenomena. For electromechanical systems, on the other hand, 
the non-relativistic form of kinetic energy must be retained for the mechanical 
degrees-of-freedom. Potential energy contributions due to elasticity or gravita- 
tion may also have to be included. 



3 Linear electric circuits 

The equations governing linear electric circuits are presented in almost every 
textbook on electromagnetism, and their similarity with those for oscillating 
mechanical systems is often pointed out. A smaller number of more advanced 
texts even go as far as presenting a Lagrangian formalism underlying the circuit 
equations [251 [^ [^. 

Here we will derive and discuss some standard results in for linear elec- 
tric circuits starting directly from ([T5|) . These are alternatively called current 
circuits, or networks, in the literature. Assume that all current flows in con- 
ducting thin (filamentary) wires and that there are n such wires with currents 
ik = Sk, (fc = 1, . . . , n). It is then easy to show that the magnetic part of (jl5p 
can be written, 

„ 1 " " 

" 2 / c^' ■ ^ " 2 ^ ^ ^^^^ 

^ k=l 1=1 

In a similar way for a fixed arrangement of m (extended) conductors, with 
charges (i = 1, . . . , m) on them, the electric part of Clc can be written, 

i=i j=i 

The inductance coefficients L^i and the generalized capacitance coefficients Fy 
only depend on the geometry of the arrangement. These are derived and ex- 
plained in the Appendices. We have thus found that the Lagrangian ([T5|) under 
the above assumptions can be written, 

CLcie,e)^ CL{e) + Cc{e). (18) 

This is a valid total Lagrangian for a non-radiating arrangement of current car- 
rying thin wires and extended charged conductors. This separation of magnetic 
and electrostatic effects comes from the central idea that there will be no net 
charge density on thin wires, and that currents in extended conductors have 
negligible magnetic effects. 
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3.1 The conductor pair condenser 

In practice all the charges ei on all the m different conductors are not indepen- 
dent. Often a circuit is arranged so that the conductors come in pairs that are 
very close, so called condensers. If each such pair is connected by a wire while 
being electrically isolated otherwise, the total charge on that subsystem must 
be a constant which we take to be zero. The number of wires is then half the 
number of conductors, n = and the charges Ci come in pairs that are equal 
and opposite ek — —et+n, {k — 1, . . . ,n — ^), while the current in the wire 
connecting them is ik = et, see Fig. [TJ There is then only n = ^ degrees-of- 
freedom of the problem. We now assume, without loss of generality, that the 
coefficients are symmetric in the indices ij, and define the new symmetric 
matrix, 




Figure 1: A single LC-circuit with an inductance L and a capacitance C. The 
Lagrangian (|20[) gives the dynamics for n interacting circuits of this type. 

Cj./ = Ci^^ =Tkl + Tk+nl+n — ^kl+n " ^Ik+n, (19) 

where k,l = 1^ . . . ,n. Using this our Lagrangian psp can be written, 

^ n n 

^Lcie,e) = 2 i^kiekei - Cj^^^ekei) . (20) 

k=l 1=1 

Here the n charges ek and currents ek are independent and a diagonal element 
of the L matrix, Lkk, is a self- inductance, while the off diagonal elements cor- 
respond to mutual inductances. A diagonal element, Cf^^ = ^kk + ^k+nk+n — 
2Tkk+n, of the C~^-matrix represents the inverse capacitance, C^fc, of the cor- 
responding conductor pair condenser. 

3.2 Equivalence of electric and mechanical oscillators 

The Lagrangian PH)) is completely equivalent to that of a mechanical system 
of coupled oscillators, the L-matrix corresponding to the mass matrix and the 
^-matrix corresponding to the stiffness matrix (of spring constants). This is 
often regarded as a purely formal correspondence, a mere mathematical mapping 
of one problem on another physically completely different one. This is wrong. 
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however. If we denote the hnear density of conducting charge in wire A; by A^, 
and the arc length along this wire by Sk, we find that the current in the wire is, 
ik = (ik = ^kSk- Here, of course, Sk is the speed of the conducting linear charge 
density. Clearly the charges on the condensers have to be, Ck ~ A^Sfe (with a 
suitable choice of origin and orientation for the arc length). If this is inserted 
in £lc we find that, 

^ n n 

^Lc{e, e) = 2 X! X! ^'^"^^ {LkiSkk - C^i'SkSi) = 
fe=i 1=1 

(21) 

^ n n 

2 5Z X! (^^ki-hsi - KkiSkSi) = T{s) - V(s) = £(s, s). 

k=l 1=1 

The generalized coordinates Sk now have dimension length so we have an or- 
dinary mechanical coupled oscillator Lagrangian on the right hand side. The 
difference is that the mass matrix, Mki = Xk^iLki, does not come from rest 
mass but entirely from the inertia contained in the energy of the magnetic field. 

By means of the technique of simultaneous diagonalization of two quadratic 
forms one can find a linear transformation to, so called, normal mode coordinates 
and thus decouple the equations of motion, which are, 

n 

J2 i^kiei + Cfe^^e,) =0, = 1, . . . , n, (22) 
1=1 

assuming that no further generalized forces enter the problem. In terms of the 
normal modes qk the equations of motion become, 

Qk + T^^fe = 0, k = l,...,n, (23) 

^k'-^k 

SO these oscillate independently with angular frequencies ujk — 1/ yLkCk- For 
the corresponding mechanical problem one finds ojk — \J Kk /Mk ■ The expres- 
sion Lo = l/y/LC for the angular frequency of a single LC-circuit, as shown in 
Fig. [1] is sometimes referred to as Thomson's formula. 

3.3 Introduction of resistance and external voltage 

Our Lagrangian (j20p corresponds to a coupled system of undamped electromag- 
netic oscillators. In most cases of practical interest the connecting wires will not 
be perfectly conducting. There will be resistance in the system. The energy will 
then dissipate and the equations of motion require that there are generalized 
forces that describe this. Assume that the ohmic resistance in wire k is Rkk- 
This can be achieved with a Rayleigh dissipation function. In a more general 
case there may also be off diagonal elements Rki , and, 

k=l k=l 

is the most general form of this function for linear circuits. 
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If there is resistance currents eu eventually dissipate to zero and arbitrary 
initial conditions only lead to transient dynamics. In most applications of circuit 
theory one is therefore mainly interested in systems with added external e.m.f.. 
This can be done by the additional term, 

n 
k=l 

added to the Lagrangian £lc- Here Vk{t) is an applied external voltage. A 
constant e.m.f. will not drive a stationary current through a condenser so to 
get a direct current some of the C~ ^-matrix eigenvalues must be zero. With 
harmonically oscillating e.m.f., Vk{t) = Vfcsin(wt), capacitors are no problem 
and one is dealing with alternating current circuits. 

The general Lagrangian equations of motion for a system of circuits are then, 

d dCcc dCcc dTZ 

:tt7Tt^ - IT^ = -7T-, 26 
at dck oek ock 

where, 

£cc(e, e, t) = £Lc(e, e) + £cmf (e, t), (27) 

is the circuit Lagrangian [24 . The system ([221) of equations of motion arc then 
modified so that, 

n 

J2{Lkiei+Rkiei + C^i^ei) =Vk{t), k^l,...,n, (28) 
1=1 

is their new form. This is thus the type of system investigated in linear circuit 
theory (see e.g. Guillemin [IS] or Josephs [M|)- In mechanical systems the ohmic 
resistance terms correspond to dampers (dashpots) and the external e.m.f. to 
applied external force. 

3.4 Energy and Hamiltonian for conservative systems 

For Lagrangians C{q, q) with no explicit time dependence, such as those of Eqs. 
([TT|) and ([2D|), the quantity, 

k=l * 

is known to be a constant of the motion, the energy. For example the Darwin 
Lagrangian (jlip corresponds to the conserved energy, 

1^/1 \ 

£d = + 2 X! ( "'^"■Wa • Aa{ra) + ea(t)a{ra) \ ■ (30) 
a=l ' 

This expression for the energy goes up if currents are parallel since the vector 
potential is proportional to terms like e^Vh- This may seem odd since we find 
in Sec. 14.31 that parallel currents attract. We will return to this in Secdjbelow. 
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Returning to circuits we find that, when there is no time dependent forcing 
and no ohmic resistance, the Lagrangian is such that there is a conserved energy. 
Using PU)) and ((^ gives the expression, 

^ n n 

Shcie, e) = 2 51 51 i^kihei + C^i^etei) , (31) 

k=l 1=1 

for this energy. RecaU that if £ = T - then £ = T +V. The effect of a 
constant forcing, due to permanent constant charge on condensers, is only to 
shift the equiUbrium from Ck — ek — Ignoring this ([20|) is the most general 
circuit Lagrangian that conserves the energy (|3ip . 

The generalized momenta obtained from the Lagrangian (j20p are by defini- 
tion, 

1=1 

The Hamiltonian is obtained by eliminating the generalized velocities in the 
Lagrangian energy (|3ip in favor of the generalized momenta. Since, Ck — 
Sr=i -^ki^Pi' Hamiltonian will depend on the inverse of the L-matrix. For a 
system of coupled LC-circuits we find the Hamiltonian, 

^ n n 

T-(-Lc{e,p) = i^ki^Pf^P' + ^ki^'^kei) , (33) 

k=l 1=1 

representing its conserved energy as a function of phase space variables. The 
Lagrangian and Hamiltonian above, as well as the interpretation of the canonical 
momenta as magnetic fluxes discussed below, can be found in an article by 
Meixner [27], discussing thermodynamic issues. 

3.5 Generalized momenta and magnetic flux 

To find the meaning of the generalized, or canonical, momenta in this case we 
return to the definition of the magnetic Lagrangian, 

£L^l[-j-AdV. (34) 



The volume integration is only over the filamentary wires that carry the currents 
ik = e/c, so, using j dV = ifcdr, we get, 

1 " r 

Ci^ = — Y^ikf A-dr, (35) 

k=l •'^ 

where the line integral is around the loop of wire k. Now, however, 

i A-dr= / (V X A) • ds = B ■ ds, (36) 

Jk Jk Jk 

according to Stokes' theorem. By definition this is the magnetic flux, <i>fe, 
through the loop k. We thus find that, 

k=l 
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Comparing with p6p this gives us that, 

(J) 

— =y^Luk=Pk, (38) 

where the individual terms on in the sum represent contributions to the flux 
through k from the loops I of the system. 

Finally then, we have found that the generalized (canonical) momenta pfc, 
of Eq. ((5^ . conjugate to the charges on the condensers are the magnetic 
fluxes through the currents loops (divided by c): pk = ^k/c. This means that if 
one e; does not appear in the Lagrangian (because there is no condenser in the 
corresponding loop) then that generalized momentum pi (or flux) is a constant 
of the motion. 



4 Electromechanical systems 

Few texts derive of equations of motion for electromechanical systems from the 
fundamental Lagrangian for particles and fields, only Neimark and Fufaev |28| 
come close. As should be clear from the above developments electromechanical 
systems, as opposed to electric circuits, require that we retreat from (|15p back to 
the Darwin Lagrangian in the form ([7|), which we had before we neglected rest 
mass inertia. From there the equations of motion for electromechanical systems 
can be found by adding constraints, or assumptions, about the motion, thereby 
reducing the number of degrees-of-freedom, in the way familiar from analytical 
mechanics. 

If macroscopic matter moves one must, of course, add the Lagrangian cor- 
responding to that motion. Further, the induction and capacitance coefficients 
may now depend on the mechanical degrees-of-freedom corresponding to the 
motion of thin wires and extended conductors of the system, since this changes 
its geometry. For an energy conserving system one then typically arrives at a 
Lagrangian of the form, 

^ n n 

jr{q, e, g, e) = T{q, q) ~ V{q) + 2 E E [Lki{q)ekei - C^i\q)ekei] , (39) 

k=l 1=1 

and this will be general enough for our purposes. One notes that if charged 
conductors move this produces magnetic effects which may have to be handled. 
In many cases, however, the speed of this motion will be such that the magnetic 
effect is negligible. In more general cases one may, of course, also have coupling 
terms between q and e. In case of doubt the safe method is to start with the 
Darwin Lagrangian (fTTj) and introduce relevant constraints and idealizations. A 
couple of examples of this procedure can be found in Essen [211 [30] ■ 

Electromechanical systems are treated e.g. in the books by Neimark and 
Fufaev [2S], Wells [H], and Gossick [32] • Articles discussing various aspects of 
these systems are [331 1311 133 [3S] • We now proceed to some concrete examples 
of magnetomechanical systems. 

4.1 Extension of current carrying ring 

When a current i flows in a conducting circular loop, or ring, its radius will 
increase somewhat. This is due to the reaction forces to the forces needed to 
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bend the current that, due to inductive inertia, otherwise would move in a 
straight Hne. Let us calculate this increase in radius. 



i>0 




Figure 2: Notation for the ring extension problem. When the current i = e flowing 
in the ring is zero its radius is R. When the current is increased the radius becomes 
R + SR = R{1 + ^). 

The self-inductance of a ring made of thin wire of circular cross-section is, 
see Appendix lA. 31 below. 

Lc{R) = 47ri? [ln(8i?/A) - (7/4)]/c2, (40) 

where R is the ring radius (at zero current), and A is the radius of the thin wire. 
Assume that the wire ring can be treated as an elastic with stiffness k. The 
energy required to increase its length, £ = 2ttR, by 6£ — 27t SR is then, 

ViSR) = hiSiy = hi2nr{6R)'. (41) 

If we introduce the notation, k — 41n(8i?/A), and, ^ = 6R/R, for the relative 
change in radius, we get the self-inductance, 

Le(e) = ^i?(l + - 7) + 4 ln(l + 0]/c^ (42) 

for a ring as a function of the relative extension ^ of the radius. The elastic 
potential energy is, 

V{0 = ^km'^\ke'e, (43) 

where k is the stiffness, or spring constant. We can now study the this two 
degree-of-freedom magnetomechanical system, the degrees-of-freedom being e, 
with i — e, and ^. 

The Lagrangian will have the form given in Eq. (|39p and becomes, 

^(e, t e) = \mR^i^ ~ y(e) + ^Lc(C)e^ (44) 

where the first term is the kinetic energy of the ring, of mass m, due to time 
dependent radius. The corresponding Hamiltonian, including the kinetic energy 
of radial ring oscillations, the elastic energy, and the magnetic energy, is thus. 

Here p = L^e is the magnetic flux, divided by c, through the ring due to the 
current e = i. The generalized momentum, = mR^^, is the momentum of 
corresponding to radial motion of the ring. 
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Assume that we increase the flux, p, through the ring slowly without exciting 
radial oscillations. We might think of it as lying on a horizontal lubricated 
surface that dissipates kinetic energy keeping = 0. Let us calculate the ring 
extension. It will correspond to the minimum of the sum of the elastic and the 
magnetic energies {£ = 2ttR), 




-OA 0.1 0.2 0.3 



Figure 3: The function E{^,p) of equation (|46|) plotted for R = lOOA, so that k « 27, 
and with, fc = ^ = p = c=l, asa function of ^. 



since this minimum corresponds to the minimum of the effective potential for 
the ^-motion, p being a constant of the motion. A graph of this function is 
shown in Fig. [31 

We first differentiate E(^,p) with respect to ^. To find the minimum we 
then wish to solve, dE/d£, = 0, for ^, but this equation does not give any simple 
analytic root. We therefore first expand in the presumably small parameter ^ 
and keep the constant and the linear term. The resulting equation is trivial to 
solve. Some algebraic rewriting make it possible to write the root in the form, 

f ( )= (^-3)(/t-7)cV 

(k- 7)3£3fc + 2[7+ (k-4)2]cV' 

where k = 41n(8i?/A). This is thus, to first order, the relative extension, ^ = 
5R/R, of a conducting ring of radius R, cross-sectional radius A, and elastic 
constant (stiffness) k, through which the current, i — p/£c, flows. The ring 
extension problem is also treated in Landau and Lifshitz, vol. 8 [24j . but in a 
more complicated way. 

4.2 Parallel coaxial circular current loops 

The mutual inductance of two parallel coaxial rings, of radius Ri and i?2, a 
distance 2 apart, is given by (Becker [57]). 

X,. . ^ r ^^^-^"^^^^ (48) 

c2 Jo ^ ^2 ^ ^2 _ 2R1R2 cos (p 
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Figure 4: Notation for tire two ring interaction problem. It is assumed that they can 
slide freely along the z-axis remaining parallel and coaxial. They have radii R\ and 
Ri and the distance between their planes is 2: = 2:1 — 22- The currents in them are 
jfc = fifc, fe = 1, 2. 



The integral can be evaluated exactly in terms the complete elliptic integrals, 
K(A:) = / , • ^ , (49) 

ri'^ , 

E(/c) = / ^Jl-k^siT^^ede. (50) 
Jo 

Putting, 



^ - V (i^TTWTI^ ^^^^ 



one obtains, 



(52) 



for the mutual inductance. 

Assuming that the rings can slide along the z-axis we now get the Lagrangian 
C = C{zi^ Z2, ii, 22, ei, 62) of this four degree-of- freedom system in the form, 

C = ^mizf + ]^m2zl + ^Lciej + ^Lc2el + ii2(zi - 22)6162- (53) 

We first do the well known transformation to center of mass and relative co- 
ordinates: Z = (mizi + m2Z2)/M,z ^ z\ — Z2, where M — m\ + TO2- The 
Lagrangian L ~ L{z, i, Z, ei, 62) is now, 

C = ^MZ^ + ^mz^ + ^Leie? + iic2e2 + ii2(2)eie2, (54) 
where m — 1111^12/ M is the reduced mass. We now study this system. 

4.3 Force per length between parallel constant currents 

First, assume that, we maintain constant current in both rings: — Ik, k — 
1,2. We note that this is a holonomic (integrable) constraint since it can be 
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integrated to give, ek — 6^(0) + Ikt, for the generalized coordinates. It is thus 
holonomic, but not time-independent {i.e. not scleronomic in traditional termi- 
nology). Since we will mostly introduce this constraint here for cyclic (ignorable) 
coordinates, i.e. coordinates that do not appear explicitly in the Lagrangian, 
the time dependence of the constraint will not be manifest in the appearance of 
Lagrangian. This is a peculiarity of the type of system treated here, which thus 
may formally appear conservative, even though it is not physically conservative. 
External energy is normally needed to maintain constant current, even for ideal 
conductors. 

With this constraint there are only two degrees-of- freedom, Z and z, and 
the Lagrangian ([54|) gives, 

C ^ ImZ^ + KnP + hhLi2iz), (55) 

where we have discarded the two constants due to the self inductances. This is 
now a simple system in which the center of mass Z-motion is trivial, and where 
Vi{z) — —Iil2Li2{z) acts as potential energy of the z-motion. Assuming that 
Ri ~ R2 — R and that z <ti R we find from expansion that. 



Ll2{Z) « 







m 


- 2 







—mz+ constant. (56) 



We throw away the constant term and find, in this approximation, the potential 
for the relative motion, 

Vi{z)^^hl2lnz. (57) 

This means that, 

£ 2nR dz c2 z ^ ' 

is the force per unit length between the rings (of length £ — 2t:R) assuming 
z <^ R. One recognizes this as the standard expression for the force per length 
between parallel currents. The minus sign means that it is attractive when 
/1/2 > 0, i.e. for parallel currents, otherwise repulsive. This force between 
parallel current carrying wires has been much discussed in the pedagogical lit- 
erature [351 SOI im HH US] but the analytical mechanical approach presented 
here does not seem to have received much attention. The importance of the 
problem originates in the fact that the definition of the ampere, the SI unit of 
electric current, is based on this type of force measurement. 



4.4 Relative oscillation of two rings of current 

Assuming constant currents is not natural in this type of problems where we 
assume energy conservation. In general maintaining constant current requires 
that work is done by an external e.m.f.. For two perfectly conducting rings of 
modest size it is more natural to assume an isolated system of constant energy. 
We return to the Lagrangian (|54p and note that since the coordinates (charges) 
ei and 62 do not appear, the corresponding generalized momenta, pi and p2, 
given by, 

dC dC 

-—=pi=Lciei+Li2{z)e2, 77^=^2 = ^0262 + ^12(2)61, (59) 
oei oe2 
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are conserved. Solving these for the currents, 



Lc2Pl- Li2{z)p2 . LclP2 - Li2{z)pi 

ei = — — ■? pr^^7 62 — 



LciLc2 — L\2{z) 



(60) 



we can proceed to find the Hamiltonian corresponding to the Lagrangian ([5^ 
This gives, 

p\ vl Lc2p\ + LcIpI - 2Li2{z)piP2 



n = 



2M 2m 



(61) 



We now wish to compare the interaction of the two rings for the case of constant 
currents ei = (=2 = /, and for the case of constant momenta, 



Pl = [Lcl + il2(0)]/, P2 = [Lc2 + il2(0)]/, 

assuming currents / at z = 0. 



(62) 




Figure 5: Plots of interaction potentials of two rings of current, see Fig. |4l For 
constant unit currents the potential is Vi — — I/i2(C) a-nd this is the upper curve. For 
constant generalized momenta, assuming unit currents at 2 = ^ = 0, the lower curve, 
corresponding to Vp of Eq. (f64|l . is obtained. 



For constant currents the interaction potential is simply the negative of the 
last term of ((5i)) . 

Vi{z) - -/2Li2(z). (63) 
For the Hamiltonian ([^T|) the potential that goes to zero at infinity is obtained 



by subtracting the constant 
potential. 



' 2Loi 



j2-, from the last term. This gives the 



2Lo 



Vpiz)^Vi{z)- 



£12(0) 

ic2 



Ll2(z) 

2 



1+ 



J-12(0) 



id 



Lc2 



LclLc2 



(64) 

for relative z-motion of the closed conservative system of two perfectly conduct- 
ing rings. 
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To get definite results and compare the two expressions we now introduce 
specific values of the parameters. We use PU)) for the self-inductances of the 
rings, Lcj, j = 1,2, and expressions (|5ip and (j52p for the mutual inductance 
Li2{z) of the two rings. For definiteness we put, i?i = £, i?2 = (9/10)^, A ~ ^/20, 
and use the dimensionless distance, = z/£, instead of the distance z between 
the planes of the rings, see Fig.[31 With these choices, and putting I = i = c = 1, 
we get the interaction potentials shown in Fig. [51 The difference is fairly small 
due two the fact that the self-inductances are an order of magnitude larger than 
the maximum value of the mutual inductance even though the parameters have 
been chosen to maximize the difference. A system for which the difference be- 
tween a (time dependent) constant current constraint and a closed conservative 
system is, qualitatively and quantitatively, of importance is described in the 
following section. 



4.5 Rectangular circuit and the rail gun 



Here we study a system which can be thought of as an idealized rail gun. It 
reinforces the lesson in the ring extension example above in showing that closed 
loops of current tend to expand. In the example of the ring extension this can be 
understood as due to the inertia of the current. The current tries to go straight 
but has to follow the conducting wire and thus there must be a reaction force 
from the current on the wire that tends to straighten it out. Here we will see 
that this also happens when there are right angled corners in the circuit. 

The self-inductance of a rectangular circuit with side lengths a and b made of 
flat conducting strips of width A in the plane of the rectangle has been calculated 
by Bueno and Assis [44l|45]. Their result is, 



Lr(a,6, A) 



1 



a arcsinh { — ] + b arcsinh 



aln(^)+61n(| 



+8v/a2 + 62-2(a + 6)|. 



(65) 



where the neglected terms are of order X {0[{X/a)^] + 0[{\/b)^]Y We now 
introduce a — x — b^ and A = b/10 and assume that ^ > 1 so that we can 
expand in 1/^. This gives. 



Lm « ^ [4 In(lO) + 6]^+ [41n(20) - 2] - - 



1 1 



e 24^3 



(66) 



where terms of order 1/^'^ and higher have been neglected. The two expressions 
5]) and (j^]) are compared in the plot of Fig. [5] 
The Lagrangian of the two degree-of-freedom system shown in Fig.[7]is then. 



C^imbH^ + ^L,{Oe^ 



(67) 



where, e = z, is the current in the circuit. The coordinate e is absent (cyclic) so 
the generalized momentum, 



(68) 
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0.2 0.4 0.6 c. 0.8 1 



Figure 6: Comparison of the two approximations of the self-inductance Lr(C) of a 
rectangle. Here h = c = 1 and the lower curve is the approximation in (I65|l . The 
upper curve is the truncated series approximation in (I66p . The approximation is 
excellent for 0.5 < i^. 



is conserved, assuming perfectly conducting parts. The conserved Hamiltonian 
becomes, 




Figure 7: The self-inductance Lr of this rectangular current loop is given by Eq. (I65p 
with X = a. K simple model of a rail gun is obtained when the right hand edge of the 
rectangle can slide, with negligible friction, in the direction of the x-axis. 



1 p1 



1 



P 



1 P| 



2m62 2Lr(0 2m62 



(69) 



where = mh^^. Assume that ^(0) = at ^(0) = so that the conserved 
energy, Ti. — E, is E 



One then finds that. 



(70) 



is the velocity of the moving bar as function of ^ = x/b. Since the self- 
inductance goes to infinity with ^ the limiting velocity of the bar will be v{oo) = 
p/ y^mL,-(^o)- The speed of the rail gun projectile is thus proportional to the 
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conserved magnetic flux $ = cp and inversely proportional to the square root 
of the self-inductance of the initial rectangle Li-{^o). 

Assume instead that we maintain a constant current, i.e. introduce the con- 
straint e = /, in the rectangular circuit. A look at the Lagrangian ([ST]) then 
shows that the there is, formally, a conserved energy, 

SitO^lnibH'-ll'UO- (71) 

The bar will move in the potential Vi{£,) = — i/^Lr(C) ^^nd thus accelerate indef- 
initely (until the end of the rails). Though this still is formally a conservative 
system it is clear that energy must be continuously feed into the system to 
achieve the continuous acceleration of the bar. In fact charge in the system 
increases linearly, e{t) = e(0) + It, according to the implied time dependent 
constraint. 

The treatment of the rail gun above is essentially new as far as the author 
knows, but this type of system has certainly been discussed both in the ped- 
agogical and technical literature (see e.g. Knoepfel [H]). Some examples from 
pedagogical journals are [47l HH |49]. Because of the many potential applica- 
tions of rail guns there is a huge technical literature on the subject. A technical 
treatment using the Lagrangian formalism of electromechanical systems is by 
Hively and Condit [5D] . 



5 On the nature of magnetic energy 

It should be clear from the above examples that the methods of analytical 
mechanics can be quite useful in treating electromechanical, and in particular 
magnetomechanical, systems. There is no need to first find the fields and then 
the forces from these. Instead both steps are integrated into a single formalism. 

While a Lagrangian with no explicit time dependence corresponds to a con- 
servative system, one should note that a constant current constraint may only 
be formally energy conserving. External work may be needed to keep current 
constant. Consider the magnetic (inductive) part of the energy expression ([3T|l . 

^ n n 

£h{q,e) ^ -^^Lki{q)ekei, (72) 

fc=i 1=1 

where q represents mechanical degrees-of-freedom. Should the currents be kept 
constant, = ~ constant, we find that this term becomes the negative of 
an effective potential for the q-motion, 

^ n n 

Mq) = - 2EE^'='(9)^fc^' = (73) 

k=l 1 = 1 

in a full Lagrangian of the form Ci{q,q) = T{q) — V/(g). For this Lagrangian 
an equilibrium position corresponds to a minimum of Vi{q). Evidently this 
corresponds to a maximum of the energy ((7^ . We have thus arrived at the result 
that in electromechanical systems, for which current is kept constant, magnetic 
energy will tend to a maximum, when the system tends to its equilibrium. 

The above result does not seem well known, even though explicitly stated in 
the textbook by Greiner [5T. It is also in accord with Woltjer's [S2] assumption 
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that self organized states of a plasma should correspond to maxima of magnetic 
energy. Mehra and De Luca |53j , on the other hand, made computer simulations 
of a plasma minimizing the velocity space form of the Darwin energy pop . They 
then found surprising results indicating that anti-parallel currents attracted each 
other. From our example above we know that it is the other way around. In 
conclusion, the velocity space form of the magnetic energy tends to become 
maximized. 

Why would magnetic energy be so different from other forms of energy which 
usually tend to minima in equilibrium states? We can resolve this conundrum 
by recalling that there is also the Hamiltonian form of magnetic energy. The 
Hamiltonian corresponding to the Darwin energy is discussed in |54j . For circuits 
we have, 

^ n n 

nL{q,p)^^Y.Y.^kl\l)PkPl- (74) 
k=l 1 = 1 

We learned that the generalized momenta pk, proportional to magnetic fluxes, 
are conserved if the corresponding charge is cyclic. For constant momenta 
Pk = Pk the effective potential for the g-motion is, 

^p('Z)-WL(g,P°). (75) 

This form of the magnetic energy is minimized when the q-motion maximizes 
the inductive coefficients. The Hamiltonian and the Lagrangian formalisms thus 
correspond to two different energy concepts. In the case of magnetic energy, 
which is a kinetic energy, the coefficients of the generalized velocities, and the 
coefficients of the generalized momenta are each others inverse, which means 
that the minimum of one is the maximum of the other. 

Schwinger et al. [6\ is the only text that discusses the above facts concerning 
magnetic energy briefly. They also point out that the velocity space form of 
energy seems to give the erroneous idea that parallel currents repel. Being clear 
about which variables are held constant when one searches for a minimum is 
essential. Stating that magnetic energy is a maximum or a minimum for some 
configuration, as has been done |55j . is meaningless unless this is made explicit. 

A Appendices 

Here the generalized capacitance coefficients and the inductance coefficients of 
linear circuit theory are derived and explained. Several results here are new. 
Even if the existence of these coefficients is stated in many texts, the discussion 
of their properties and physical origin usually is rather brief. 

A.l Energy of a system of charged conductors and gener- 
alized capacitance coefficients 

Consider the electric energy expressed in the form, 

W,^^Jg<j)dV. (76) 

The potential, (p, is a solution of Poisson's equation, 

= -Attq, (77) 
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where it is assumed that the solution must go to zero at large distance from the 
region where the charge density, g, is located. For point particles the solution 
to this equation is well known. Assuming that, q = Qi(r) = ei6{r — Ti), the 
solution is, (/) = (t>i{r) ~ ei/\r — r^j. For N particles the charge density is the 
sum, Q = gi, and the solution is simply a superposition of such solutions, 
= X^f^'/'jj "i^s to the linearity of the Poisson's equation. When these 
results are inserted into (|7^ one obtains the result. 



where, — \rj — r^j, are the distances between the particles, and where it is 
necessary to exclude the case, i = j, to get a finite result. 

Assume now that we, instead of particles, have a set of Uc fixed conductors 
occupying the (compact) volumes 14, (fc = l,...,nc) and ask: what is the 
potential (f> and the energy We of the system if we put charges e^- on all, or some, 
of these (isolated) conductors? We know that the charge on each conductor must 
be distributed on its surface in such a way that the electric potential, (l){r), is 
constant on the surface, and the interior, of each conductor 14 . Otherwise 
current will flow until it becomes constant, since a gradient implies presence of 
electric field. 



Figure 8: Here we study the potential arising from a system of conductors Vi which 
either have a distribution g*^ of unit net charge or are only polarized with zero net 
charge redistributed to a density Q*j, so as to make the potential constant on each 
conductor. In the figure i,j = 1, 2, 3. 

For the particle problem the charge density, g = Sig* , with, g* = S{r — 
Ti), gave rise to the potential, = ^i^l^ where, </>* = l/|r — r^j. We now 
wish to determine corresponding densities g* and potentials (j)* for the case of 
isolated conductors, of given size, shape, and position. Assume that we place 
a unit charge on conductor Vi while all the other conductors have zero charge. 
We first note that this does not mean that the charge densities on the other 
conductors are identically zero since charge on one of the conductors will induce 
a polarizing redistribution of the charge on the others. The total charge density 
for this case is thus. 



1 



N 




(78) 






(79) 
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where, 

9:u^V-5.,^[l I 1 = 1 (80) 

i.e. the net charge is unity on conductor i and zero on the other conductors. 
Our assumption means that these densities must obey, 

V^c/,* = -Ang*, (81) 

with (j)*{r) given by, 

and being such that, 

(/)*(rj) — Tij = const, for aU rj E Vj. (83) 

Properly chosen g*j must thus produce functions (/)* are constant on all of the 
conductors. One notes that the field (puir) has monopole character, while the 
(f>*j{r), for i 7^ j, are of a dipole character. 

The general solution to the Poisson equation ([77]) for a system of Uc conduc- 
tors with charges on the first A^(< ric) of them is then, 

N 

c^{r)=J2e,cl,:{r), (84) 
1=1 

according to the superposition principle. We see that the potentials 4>*, arising 
from unit charge on conductor i and zero on the rest, constitute a basis set of 
functions for this problem, analogous to the functions 0* = l/|r — r^j for the 
point particle problem. 

Let us return, now to the energy expression Eq. (j76p . Since the charge 
density g is non-zero only on the conductors Vk we find that 



l f Q^dV^^Y. f g{r)^{r)dV. (85) 



Use of (|84| gives us, 



1 r 



but, since 0* is constant on each of the conductors it can be taken out of the 
integral and replaced by its constant value on conductor Vk, 

, Uc W „ N N 

aEE^^'^Hr-;.) / e(r)dF=-5]5]e.</.r(r,)e,. (87) 

fe=l i=l ''^k j^i i=i 

The second equality here is due to the fact that, 

N N 

^? = E ■^'j = E E 

3=1 i=i 1=1 



21 



so, according to ([50]) . the integral of the charge density over Vk is simply Ck, for 
k — 1, . . . ,N, and zero for fc = + 1, ric. So, using we finally get, 

1 ^ 

ij = l 

for the energy of the set of conducting bodies. When we compare this to ([78|) we 
see that the quantity represents a kind of effective inverse distances between 
the charges on the conductors. As defined in ([5^ - [55)1 this quantity is the 
constant value of the potential on body j when unit net charge is distributed 
on body i and zero net charge on all the other Uc — 1 bodies. Note that bodies 
with zero charge contribute to the values of the Tij. 

If we denote the constant value of the potential on conductor i by (pi — 4>{ri) 
we find from ([M]) that, 

N N 

(|y^=Y. = ^1^3^^ for i = 1, . . . , iV, (90) 

i=i i=i 

where ([55]) was used to get the last equality. These linear equations for the 
potentials can be solved for the charges. The result can be written, 

N 

e, = "Yaj4)j, for i = l,...,iV, (91) 
i=i 

where the coefficients Cy represent the matrix elements of the inverse of the ma- 
trix with elements Tij. The Cij are called generalized capacitance coefficients. 



A. 2 Energy of a system of current carrying wires and the 
inductance coefficients 

This is a subject discussed in many textbooks. Examples are Landau and Lif- 
shitz, vol. 8 [23, Greiner's Classical Electrodynamics and Johnk ^Bj. A 
more specialized text is by Knoepfel [IS] . 
The magnetic energy is given by, 

W^^^^jj-AdV. (92) 

If one assumes that all current is fiowing in thin wires one can replace the 
infinitesimal vector j AV with dr^ where dr^ is a line element along the 
curve Cfc defined by the fcth wire (filament). This gives. 



k 



One notes immediately that if the curve is a closed curve the corresponding 
contribution is gauge independent since, 

/ A-Aru^i (A' + Vx) • drfc = / A' • dr^. (94) 
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For conducting wires that are not closed but instead go between conductors that 
provide capacitance the issue of gauge may need to be resolved. 
When there is gauge independence we can use the expression, 

for the vector potential from the currents ii flowing in wires along the curves 
Ci. Inserting this into we find that the magnetic energy is, 



if Jcjc, \rk-ri\ 
So, if we introduce the wire geometry dependent quantities. 



Jcjci \rk-ri\ 



we find the magnetic energy in the form, 

t^m = i^ifciJfeii. (98) 
kl 

For fixed positions of the wires the magnetic energy is thus a quadratic form in 
the currents with constant coefhcients Lki- For I ^ k these are called mutual 
inductances. When / — k they are self-inductances. 

Apart from the problem of lack of gauge invariance for non-closed wires one 
must also deal somehow with the logarithmic divergence of the self-inductance 
for a truly filamentary wire. If the wire is assumed to be a mathematical curve 
of no thickness the expression (|97p will diverge when k — I. Should one find a 
way to handle the divergence there remains the question of the gauge invariance 
of the self- inductances. Let us consider these questions by means of an example. 

We first note that if the Darwin expression ^ for the vector potential, is 
used instead of (f95|) . the induction coefficients become. 

Iff dvk ■ dri + {dvk ■ eki)(dri ■ eu) 
^ ^ ic. ic, ' ^''^ 



i.e. different from the Neumann form (j97p . 



A. 3 Self-inductance of a rotating polygon of charged par- 
ticles 

One way of handling the divergence in the filamentary self-induction is to con- 
sider the current as due to many charged particles travelling in the wire, instead 
of a continuous distribution of charge. Assume that the current in a circular wire 
is due to N particles of charge e forming a regular polygon that rotates rigidly. 
We assume as usual that the material of the wire is of the opposite charge and 
cancels the charge of the particles so that effects of Coulomb interactions are 
negligible. Assuming that the circle has radius R, positions and velocities are 
given by, 

ri{Lpi) = Rep{(pi), and Vi{(pi, (p^) = Ripie^{ip^), (100) 
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Figure 9: An example of a polygonal circular current with A'" = 6. 

where, ep{ip) = cos (p + sin (p By and Cp = (pe^, as usual. If we further assume 
that all the angular velocities ipi are the same, and given by ipi = u, and that 
the angles are ipi = {i — 1)2tt/N for i — 1, . . . ,N, we have a rigidly rotating 
regular polygon. This is illustrated for = 6 in Fig. [5] 

We now calculate the magnetic energy of this system and compare the ex- 
pressions obtained using the two different expressions ^ and fSS]) for the vector 
potential. Simple calculations show that the magnetic energy can be written, 

N 

W„^ = ^yv,A^, (101) 



2c^ 

2 — 1 



where all Vi = Ripi — Rlu, and where, 

N 



Here V,^ is given by. 



^^^-^T.^iV^i'P^-Vj)- (102) 



in the Darwin case (fT5|) and by, 

v^H.) = (104) 

\/2[l — cosip) 

in the Lorenz case (|95p . Because of the symmetry all terms in the sum (|10ip 
are equal and we find, 

^'""2^(D N'i — J2^^iiA^)A^\, (105) 

where Aip — 2n/N. We note that the current in the ring is i = so we can 
write the above expression in the form Wm — \Li^ where L, by definition, is 
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the self-inductance, 



(106) 



of the polygon. 

For all but the smallest values of A'' it is useful to approximate the sum with 
an integral. Using the trapezoidal rule [57j . 



fix) dx 



b-a \ 1 



fib) 



(107) 



we find that (taking n = N — 2), 



N-l 



V^iiAip)Aip 



1=1 



27r 



27r 



(108) 



Using this and expanding the result in powers of 1/N one finds the results, 

in 









- 1 







for the Darwin case p03p . and, 



An 



-R 



In 



v) 



(109) 



(110) 



for the Lorenz case ()104p . The integral can be done analytically. Terms in 
l/N'^ and higher were neglected in the expansions. Maple [58l was used in these 
calculations, as well as for most other calculations and plots of this review. 
These results may be compared with the traditional result (Becker |37)), 



Lc — -irR 



In 



(111) 



for the self- inductance of a circular loop conductor used above in (|40p. Since 
N — 2nR/5, where S is the distance between neighboring charges along the 
circle, the logarithmic parts of the results agree ii S — A/2. Equivalently the 
polygon must have N = 4nR/X electrons to agree with the logarithmic part of 
the ring result, A being the cross sectional radius of the ring. 

For large N the logarithmic part of the self-inductance should be the dom- 
inating one. The contributions linear in R are seen to be different in all cases 
but do seem to have an order of magnitude agreement. 

There is a considerable literature on inductance calculations. The books by 
Bueno and Assis [IJ and by Grover may be mentioned. Technical papers 
are [BUI [HIl El] , and a few applications of ideas concerning self-inductance can 
be found in [Ml [Ml iS]. 
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